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^ Abstract 



We consider the energy modeling a two component Bose-Einstein condensate in 
the limit of strong coupling and strong segregation. We prove the T-convergence to 
a perimeter minimization problem, with a weight given by the density of the con- 
densate. In the case of equal mass for the two components, this leads to symmetry 
breaking for the ground state. The proof relies on a new formulation of the problem 
in terms of the total density and spin functions, which turns the energy into the 

\jZ. sum of two weighted Cahn-Hilliard energies. Then, we use techniques coming from 

geometric measure theory to construct upper and lower bounds. In particular, we 

q make use of the slicing technique introduced in [6] . 

m 
> 1 Introduction 

•i-H 

X 

The aim of this paper is to prove a T-convergence result for a functional modeling a 
two component Bose-Einstein condensate in the case of segregation. We introduce a 
new formulation of the problem which transforms the two wave functions describing each 
component of the condensate into total density and spin functions. The new functional 
in the density and spin variables is given by the sum of two weighted Cahn-Hilliard 
energies modeling phase transition problems as in the Modica-Mortola problem [26]. In 
fact, our new functional is strongly related to that of Ambrosio-Tortorelli approaching 
the Mumford-Shah image segmentation functional [6]. We use techniques coming from 
geometric measure theory [3, 4, 6, 10] to construct upper and lower bounds for our initial 
functional and prove T-convergence to a perimeter minimization problem, with a weight 
given by the density of the condensate. There is a large mathematical literature about the 
segregation patterns for two component Bose Einstein condensates [8, 9, 13, 14, 28, 30]: 



regularity of the limiting functions, regularity of the interface, asymptotic behaviour near 
the interface. All these papers use the limiting equations and do not take into account 
the trapping potentials and the T convergence of the energy as we do. 
Before introducing the functional for a two component Bose Einstein condensate, we recall 
some properties of a single Bose Einstein condensate (BEC). A single BEC is described 
by the wave function r) minimizing the energy 

EM = \ j^ |Vrf + ±V(x)\r,\ 2 + ^H 4 (1-1) 

where V is the trapping potential, usually taken to be harmonic, that is V(x) = \x\ 2 , e 
is a small parameter giving rise to a large coupling constant describing the repulsive self 
interaction of the condensate. The minimization is performed under the mass constraint 
J R 2 \v\ 2 = 1- We define the ground state by 

E £ {r le ) = p inf E e { V ) , (1.2) 

Jr2 M 2=1 

which is, up to multiplication by a constant, a real positive function. Let 



p(x) = max(A — \x\ , 0) with A > chosen such that / p — 1 where T> = B(0, A). 

Jv 

(1-3) 
Then, when e is small, the ground state i] £ is close to the function ^fp in V, with expo- 
nential decay at infinity. Properties of i] £ can be found in [1, 2, 17, 20, 21]. 
A two component Bose Einstein condensate can be experimentally realized as 2 isotopes 
of the same atom in different spin states [19] or isotopes of different atoms [27]. They are 
described by two wave functions u\ and u 2 , respectively representing components 1 and 
2. The Gross Pitaevskii energy of the two component condensate is given by 

S £ (u l7 u 2 ) = E £ ( Ul ) + E £ (u 2 ) + -g £ M 2 |w 2 | 2 , (1.4) 

where E £ is given by (1.1) and g £ is the intercomponent coupling strength. The energy is 
minimized under the mass constraints 

\uj\ 2 = oij with cij > and a± + a 2 — 1 . (1.5) 

In [24], numerical simulations have been performed to classify the ground states according 
to the values of e, g £ and also the rotational velocity. For e small and g £ large, the 
numerical evidence is that, for «i = a 2 = 1/2, the preferred ground state is such that 
each component is asymptotically located in a half disk with a local inverted parabola 
profile. If «i ^ «2, they occupy sections in a disk, the area of which is proportional to 
«j. In particular, when neither ai is too small, this configuration has less energy than a 
disk vs annulus configuration, which also provides segregation but preserves symmetry. 
Observation of symmetry breaking has also been obtained experimentally very recently 
[25]. The breaking of symmetry has been analyzed in [29] in a different limit, namely in 
the case e large and g £ large. 

Here, we assume strong coupling between components, that is, g £ — > oo, and we study 
the regime 



g £ e 2 — y +00 and e — >■ . (1.6) 

A trick introduced in [24] is to use a spin formulation also called the nonlinear sigma 
model. In our special setting, since the ground states are non vanishing real functions, 
this amounts to defining 



v:= vNl 2 + l^l 2 and ^ : = Arg ( l^l+jj ^L), (L7) 

where f] £ is defined in (1.2). The definition of ip implies that |wi| 2 — |-U2| 2 — i] 2 v 2 costp. 
The mass constraints (1.5) can be written as 

r/ 2 v 2 = oci + «2 = 1 and / r] 2 v 2 cos if = a,\ — a 2 • (1-8) 



We point out that cos if corresponds to the third component of the spin function. Because 

there is no rotation in the system, the ground states are, up to multiplication by a complex 

number of modulus one, positive functions. Thus, the second component of the spin is 

zero and the first one is sin if. 

Since the components are expected to segregate, the expected behaviour is thus that v 

tends to 1 except on a transition line corresponding to the interface between the two 

components, while if tends to on component 1 and it on component 2. This is what we 

want to analyze rigorously. 

We split the energy into its main contributions and will prove that 

£ e ( Ul , u 2 ) = E £ ( V£ ) + F £ (v) + G £ (v, ip) (1.9) 

where E £ is given by (1.1), i] £ is the ground state of E £ and 

Fe(v) = If v'lVvf + ^til-v 2 } 2 , (1.10) 

G £ (v,f) = I f rjy\Vif\ 2 + V yg £ {^-^s 2 (<f)} (1.11) 

and g £ = g £ ( 1 ^2 ) • Since r] 2 converges to p given by (1.3) in T>, the limits of F £ and 

G £ can be analyzed as the limits of 

\j v P\^\ 2 + ^ 2 p 2 {l-v 2 £ } 2 (1.12) 
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pv 2 £ \Wf £ \ 2 + p 2 vt~g £ {\ - cos 2 (if £ )}. (1.13) 



These two energies are of Modica Mortola types with a weight which vanishes on the 
boundary of T>. Given the definition of ip £ , there is a domain where cosip £ tends to 1 
(asymptotic region of component 1) and a domain where cosip £ tends to —1 (asymptotic 
region of component 2), and thus a transition region exists between the two domains. 
Two options exist for v £ : 

• either v £ goes to 1 everywhere, which makes the first energy small and the second 
energy of order \fg £ , 



• or v £ goes to zero on the transition line where cosy? varies from +1 to — 1: this 
makes the second energy of lower order and the first energy of order C/e. 

Because of our hypothesis that e 2 g e tends to infinity, it is the second scenario which costs 
less energy. Though v £ goes to 1 on each component, it has a transition region of size e 
where it goes sharply to zero. The second energy is of lower order and cannot be seen in 
the limit. It has just the effect of creating a small region around the interface where v £ 
is small. The first energy can be analyzed with techniques coming from [6] and, once the 
rescaling in e is made, the r-limit comes from the problem on lines: 

I(x) = inf J - / p(x)(w') 2 + -p{xf{l - w 2 ) 2 ; w G Lip(R+) ,iu(0) = and w(+oo) = 1 i 

Using the Euler-Lagrange equation associated with J, we shall see that for x G T>, the 
infimum is attained by the function 

w x (t) = tanh 
and we shall have 

1 r 1 

I(x) = ap{xf' 2 with a = — \ {l-t 2 }dt. (1.14) 

v 2 Jo 

This means that w x is the optimal profile transition at the point x, and that ap(x) 3 / 2 
is the minimum energy needed by w, to go from to 1 at x. In the ID direction, this 
provides a weight 2ap{xfl 2 because as e — > 0, v £ goes from 1 to on one side of the 
interface between the two components, and from to 1 on the other side. Therefore, we 
expect the limit to be defined as the integral on the interface where p goes from to it of 
the function 2ap[xf^ 2 ■ This requires a precise mathematical definition for this interface. 
We define X as the space of functions (p G BV\ oc (V ; {0, n}) such that 




pcos cp = «i — a 2 ■ (1.15) 

We will prove the T-convergence of s(£ e (-, •) — E e {j] e )) to J 7 given in X by 

Hv) = — [ P h |£VI • 



The limiting energy T measures the length, with a weight of p 3 / 2 , of the interface between 
the two phases of (p. Each phase of •p corresponds to one component of the totally 
segregated two-component limiting condensate. Notice that when T{tp) is finite, {p = n} 
has finite perimeter in compact subsets of T>, and 

F(p) = 2a I p^dH 1 = 2a f p^dU 1 . 

Here d*{p = tt} stands for the reduced boundary of {p = n} and Sp is the complement 
of the Lebesgue points of p, that is, 

Sp= \x eV;$teR such that lim — / \p(y) -t\dy = o\ . 

I ^o+ Trr 2 J Br(x) J 



We refer to [5, 15, 18] for the geometric measure theory concepts. We also refer to [3] for 
an introduction to the theory of T-convergence and to the Modica-Mortola theorem by 
G. Alberti. 
We now state our main theorem: 

Theorem 1.1. Let us assume that V(x) = \x\ 2 , and let 
U= \(u 1 ,u 2 ) G H l (R 2 ;R) x H\R 2 ;R) , / V(u\ + u\) < oo , («i,u 2 ) satisfies (1.5) 



The functional e(S £ (- , -) — E £ (r] £ )) V -converges with respect to the L] oc (T>) xLl oc (V) distance 

to J-{tp), in the following sense: 

(Compactness) for every sequence {(wi j£ ,W2,e)}e>o of minimizers of 8 £ in % such that 

sups(£ £ (u 1:£ ,u 2 , £ ) -E e {r} e )) < +oo, (1.16) 

£>0 

there exists tp G X and a (not relabeled) subsequence such that 



(«i, e ,«2, e ) -»• yfp (l {v ,=o}, !{*=*}) in L) oc (V) x L] oc {V) ; (1.17) 

and (Lower bound inequality) 

limmie{S £ {u h£ ,u 2y£ )-E £ { V£ )) > ?{<p) . (1.18) 

(Upper bound inequality) For every tp E X, there exists a sequence {(wi, e , W2,e)}e>o C 
"H, converging as e — > to y 7 /) (l{ v =o}) l{^=7r}) iw L} 0C {V) x L} 0C {V), such that 

limsupe(£' c (ui )e ,tt2 )e ) -E £ (r] £ )) < F(tp) . (1.19) 

We point out that we only prove the T-convergence at the level of minimizers of S £ . 
Indeed, minimizers of the functional have the property that they are positive functions 
which do not vanish. Therefore, this property allows the definition of (v, tp) through (1.7). 
As usual, the T-convergence theorem implies the convergence of the energy of the ground 
states: 

Corollary 1.2. If {(ui je ,U2, £ )} £ >o is o sequence of minimizer of S £ in %, then 

\ime(£ £ (u 1:£ ,u 2 , £ )-E £ (r} £ ))=iniF. (1.20) 

A study of the ground states of J 7 allows us to prove symmetry breaking when neither ctj 
is too small: 

Corollary 1.3. There exists 5q of order 0.15, such that if ot,\ G [Sq, 1 — So], then for e 
sufficiently small, the minimizers (ui y£ ,U2 )£ ) of £ £ in % are not radial. 

Remark 1.4. Our main theorem remains true when V is any trapping potential for which 
we have good estimates for the ground state r\ £ , namely the estimates in Proposition 2.1. 



1.1 Links with related problems 

The segregation behaviour in two component condensates has been widely studied: reg- 
ularity of the wave function [14, 28, 30], regularity of the interface [13], asymptotic be- 
haviour near the interface [8, 9]. The main difference with these references is that, on the 
one hand, we use mainly the energy instead of the equation and, on the other hand, we do 
not switch off the trapping potential by blowing up the problem near the interface or by 
considering a bounded domain with no trapping. Indeed, we consider the limit where e 
goes to zero at the same time as g £ e 2 going to infinity, so that it is the trapping potential 
which provides the leading order behaviour of the wave function through the inverted 
parabola profile p. In all the previous quoted references, e is set to 1, so that in the limit 
g £ large, the trapping potential is not present, and the limiting profile is 1. We deal with 
the trapping potential by a proper division of the limiting wave function which allows to 
express nicely the energy using a trick introduced by [22] . Nevertheless, our proofs which 
rely on energy considerations also provide information for the case p — 1. 
In [31], the authors fix a point x^ on the interface dA, and consider a sequence x £ tending 
to x^ such that u\ y£ (x £ ) = U2, £ {x £ ) = m £ . An open question in [31] is to prove in 2D that 
g £ m\ stays bounded. This may be obtained with our technique since in our case m £ is 
probably related to the minimum of v £ . We detail this remark in Section 5.3. 

1.2 Main ideas in the proof 

Let us now give more details on the proof. 

The proof consists of upper and lower bounds, that we construct for the functional 

Fe(v £ , Ve) = £ (£ e (Ml, e , U 2 , £ ) ~ E £ (t] £ ) ). 

For the upper bound, we choose the set A where asymptotically u 2 will be p. In a first 
step, we assume that (p = ttIa, where A is an open bounded subset of IR 2 with smooth 
boundary such that l-l l (dA fl dT>) = 0. The test function (p £ is matched between in a 
subdomain of T> \ A to n in a subdomain of A, using a transition region of size et £ . In 
order to approximate the optimal 1 dimensional profile that solves I{y), we define 

m £ in (0,t £ ) 

tanh in (t £ ,T) 

< r - r= \ h in (T,T + i/t) 

1 in (T + i/t,+oo), 

where t £ = tanhm e and h is a polynomial which matches smoothly tanh to 1. Then we 
define 

<At) = w £ , T (^ft 

for t = d(x)/e < CT , and d(x) is the distance to the boundary. In order to construct 
v £ , we need a partition of unity for dA, where we match the functions w v ^ T , as y,i varies 
along this partition. For this v £ , we can estimate F £ with techniques similar to those of 
Modica Mortola [26] , and to the adaptation of these techniques to problems with weight 
by Bouchitte [10]. Because p vanishes, we cannot use directly the results of Bouchitte 
and we need precise estimates on the behaviour of r] £ near the boundary. Since w £) t is 
the optimal profile for the ID version of (1.12), there is a transition from 1 to and a 
transition from to 1 and we find an upper bound which is 2 J dA I(y) dy. Then we prove 
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that for this test function, G £ (v £ ,ip £ ) is lower order: indeed, the transition layer for tp is 
is of order et £ , so much smaller than the one of v £ . Hence in G £ , v £ can be approximated 
by m £ . We choose m\ = e 2 g £ , which tends to 0, and makes G £ of lower order. 
This provides the upper bound for an open bounded subset A with smooth boundary 
such that l-i l {dA fl dT>) = 0. We show in the appendix that for any <p e X, {(p = it} 
can be approximated by sets A which are open bounded subsets of IR 2 with smooth 
boundary such that l-l l (dA fl &D) = and that the mass constraints can be satisfied for 
the approximating Ui t£ ,u 2 , £ . 

The difficulty in the lower bound is to prove that v £ goes to zero on a line and that it 
provides a positive lower bound. Indeed, the usual Modica-Mortola bound would imply 
that v £ goes to 1 almost everywhere and the lower bound is 0. We have to use G £ and the 
upper bound to prove that v £ has a transition to and that cos 2 ip e tends to 1. Hence, 
because of the mass constraint, we get two regions where asymptotically tp e is and n. 
To analyze the behaviour of v £ , we use the slicing method introduced in [6] (see also [11]). 
This consists in looking at the transition for v £ in one dimensional slices and get the ID 
energy estimate. The use of the energy G £ is only to prove that v £ goes to zero. We first 
prove the lower bound for sF £ in ID using the coarea formula, and then in 2D using the 
slicing method. We get that eF £ (v £ ,(p £ ,E) converges to a measure fJ,(E) supported in S v 
of density p 3 / 2 with respect to the "H 1 measure. The last part of the proof of the lower 
bound is inspired by ideas in [4]. 

We end with a variant of the coarea formula that can be found in [23] Lemma 2.2, and in 
[10] Proposition 2. 

Proposition 1.5. Let Q be an open bounded subset ofW N , and ^(x, s,p) a B or el junction 
o/DxKx IR^, which is sublinear in p. Let u be a Lipschitz continuous function on Q 
and denote, for every t > 0, S t = {x G fl ; u(x) < t}. Then, for almost every t 6 1 , lg t 
belongs to BV(fl) and we have 



Jn 



ty(x,u,Du)dx = / dt / V(x,t,Dlo t ). (1.21) 



The paper is organized as follows: in Section 2, we present the properties of r] £ . Then 
in Section 3, we prove the decoupling of energy (1.9) and how to go from the (wi, - ^) 
formulation to (v , ip). Section 4 is devoted to the upper bound, and Section 5 to the lower 
bound. Finally, in Section 6, we prove our main theorem. 

1.3 To go further 

1.3.1 Analysis of the limiting problem 

A natural question is to analyze the limiting problem, that is the ground state of J 7 under 
the constraint (1.15). If we define A to be the set where cosy? = 1. Then j A p = a\ and 

Iv\A P = a 2 wrtl1 «1 + a 2 = 1- 

If p = 1, then the problem of minimizing J 7 amounts to minimizing \dA\ under the 
constraints \A\ = a\ and \T> \ A\ — a 2 — 1 — ol\. The Euler-Lagrange equation of the 
minimization problem yields that the curvature is either or constant, hence A is either 
a disk, an annulus or a disk sector. The equivalent problem with a weight p is open. 
If we assume that the solution is either two disks sectors or a disk and an annulus, we can 
compute explicitly the energy J 7 and find that if a\ = a 2 , then the optimal configuration 



is two half disks, while if ot\ is much less then «2 5 then the ground state is a disk and 
an annulus (see Section 6.4). Indeed, the energy of two disk sectors is 3<r/2, while the 
energy of a disk and annulus is 8cr(l — a 1 ) 3 / 4 (l — ^1 — Oi) 1 / 2 if a.\ corresponds to the mass 
of the inside disk. If «! or « 2 = 1 - oi is to small, then the disk and annulus becomes 
the preferred configuration. In the case 0:1 = 0:2 = 1/2, it follows from our theorem that 
symmetry breaking occurs since at the limit, the disk plus annulus configuration does not 
minimize the energy. These two cases are well illustrated in the experimental observations 
of [25], figure 4. 

We insist on the point that a rigorous analysis of the ground states of J 7 in X is an 
interesting open question. 

1.3.2 Convergence for w l£ , w 2e 

The convergence that we have for (tii, e , -u 2 ,e) to y/p(l^ =0 y,l^ =7T y) is very weak. Never- 
theless, we expect that on compact subsets of l{^=7r} or 1{^ =0 }, the convergence can be 
improved. For instance, it would be natural to have similar convergence as that of i] £ to 
^fp (that is Cl oc ) on these domains. 

1.3.3 Case g £ e 2 of order 1 

An interesting open question is to deal with the case when g £ e 2 tends to a positive finite 
constant Cg. In this case, F £ and G £ become of the same order and we expect that 
m = lim inf £ _>o v £ is a positive constant (on the interface where tp varies), instead of being 
0. We believe that our techniques still provide an upper bound for the problem. We 
expect the T limit to be 

(2a m + c ^m 3 ) i f p :i ' 2 \D<p\ . 
where a m = ^= f m (l - t 2 ) dt. 

1.3.4 Case of different scattering lengths 

In this paper, we consider that the scattering lengths are the same for both components, 
that is, in (1.4) it is the same energy E e for both components. When the two components 
result experimentally from different atoms, the two scattering lengths are very close but 
not equal. This leads to an energy E £ depending on the component, namely 



E *M = \ I |v?7| 2 + 4m 2 M 



2 1 ft 1 1 
2e 2Ul 



where g^ is related to the scattering lentght of component %. li gi ^ g 2 , then the leading 
order Thomas Fermi approximation is no longer the same for each component, namely it 
is 

giPi = X 2 - \x\ 2 in B { = B(0, A*). 

The limiting problem becomes: find a partition of B\ U B2 into three sets A±, A2 and N, 
such that u 2 £ — > pi^-An J A . Pi = ®i and it minimizes 

\ x \ Pi + ttPi + / Fl P2 + -wP 2 - (I- 22 ) 



This problem is open and is probably related to the problem of finding a partition of the 

disk into two subdomains which minimize the sum of the first eigenvalues of the Dirichlet 

laplacian. 

Of course, in our case, since we have Bi = B 2 , p\ = p 2 an d N = 0, (1.22) does not 

provide any information at leading order. This is why we have to go to the next order 

which yields the perimeter minimization problem. 

2 Estimates for r\ e 

Let r] £ be the ground state defined by (1.2). The ground state is a non vanishing radially 
symmetric function. It is unique up to multiplication by a constant of modulus one, and 
satisfies the Gross-Pitaevskii equation 



- Ar] £ + —\x\ 2 r] £ + —\r] £ \ 2 r] £ = ^ r] £ . (2.1) 

e z e z e z 

The term e~ 2 X £ is the Lagrange multiplier associated with the mass constraint, and the 

pair (i] £ , \ £ ) is unique among positive solutions of (2.1). As e tends to 0, r\ £ tends to yfp 

given by (1.3). Throughout the paper, we will need precise estimates for this convergence. 

The following proposition, based on previous results in [2, 16, 17, 20, 21], sums up the 

properties of i] £ . We point out that it follows from [16, 17, 21] that an approximation 

of i] £ by yfp holds as close to the boundary of T> as needed and is given by (2.5). We 

also include an estimate of p in terms of the distance to the bulk that will be used in the 

proofs. 

Proposition 2.1. There are constants c, C > 0, a G ( 1 /2, 3 /5) and 7 G ( 1 /2, 3 /4), such that 
for e sufficiently small, p,X being given by (1.3), 

E £ ( V£ ) < C/e 2 , (2.2) 

|A £ -A| < Ce|lne| 1/2 , (2.3) 

\\Vs ~ VpWchk) < C K e 2 \\ne\ for K CC V , (2.4) 

\rj £ {x)- y/p{x)\ < Ce 1 forxeB(0,\-ce a ), (2.5) 

r] £ (x) < c e 1/G e C£ ~ 1/3{x ~ lxl) forxeM 2 \D, (2.6) 

d r 7] £ (\x\) < forxeM 2 , (2.7) 

p(x) 



G [1,2) forxeV. (2.1 



Xdist(x,dV) 
Proof: for the proof of (2.2), one can rewrite the energy as 



Ee(v)=E 1 M + ^(\ 2 - 1 - j^p 2 ) (2.9) 



2e 2 V" 2 



where 



eKv) = \l |v?f + -L (M 2 - P (x)) 2 + i(A 2 - \x\ 2 y 

1 j R 2 is e 



\v\ 2 



and (A 2 — |ic| 2 )- is the negative part of (A 2 — \x\ 2 ). In Theorem 2.1 of [2], it is proved that 
E l £ (n) < C\ lne|. Then (2.2) follows from (2.9) and the fact that J v p 2 = 2A 2 /3. 
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Estimate (2.4) is proved in Proposition 2.2 of [20]. Estimates (2.3) and (2.6) are proved 
in Theorem 2.1 of [2]. Estimate (2.7) is also proved in Theorem 2.1 of [2], but only in 
a neighborhood of dV>. But the proof, however, works in the case V(x) = \x\ 2 and the 
estimate holds in all 1R 2 . 

We now prove (2.5). For A > 0, we define rj £y \ as the unique radially symmetric, positive 
solution of the equation 

-e 2 Ar] + (\ 2 - \x\ 2 )r] + r] 3 = . (2.10) 

The function fj £i \ corresponds to a ground state of a BEC without mass constraint. In 
[16, 17, 20], the behavior of fj £y \ is studied. Using the results in Proposition 1.2, Remark 
1.3 and Proposition 1.4 in [16], we obtain 

V £ , 1 (x)=e^v (^f^+O(e), 
where 

My) = y 1/2 - \y' 5/2 + eWo(*r 11/2 ) , v e (-00, e" 2/3 ] . 

Hence, for x G -8(0, 1) we obtain 

| fj £:1 (x) - Vl-M 2 1 < C (e 2 (l - \x\ 2 )- 5/2 + e\l - \x\ 2 )- 11/2 + e) . 
In particular, if x G B(0 , A — e a ) with a G (}/2, 3 /s), we get 

|77 e ,i(a;) -y/l-\x\ 2 \ < C (e 2 e- 5a/2 + e 4 e- lla/2 + e) = 0{e<) (2.11) 

with 7 G ( 1 /2, 3 /4). We will use (2.11) to prove (2.5). First, a straight computation shows 
that defining E\ = X~ 2 e, fj £Xy \ solves equation (2.10) with A = 1. Hence, considering 
(2.11), a change of variables gives 

IUx)-Vp(x)| = 0(e 7 ), (2.12) 

for x G -8(0 , A — (\~ 2 s) a ). In Proposition 2.2 and Theorem 2.2 in [20], it is proved that 

||Vr ?£ || L o 0(M2) = 0( £ - 1 ); (2.13) 

and that 

^)=^MC>). ( 2 - 14 ) 

where 

4,A=(l + ^) and ~e = t e \e. 
It follows from (2.3) that 

4,A = 1 + C(e 2 |lne| 1 / 2 ) and e = e + 0{e 2 \lne\ 1/2 ) . 



Hence, using (2.13) and (2.14), we obtain 
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Ve,x(x) = Ve,x^) + 0(e\lne\ 1 / 2 ). 
Putting this last estimate in (2.12), and using that 7 G ( 1 /2, 3 /4), we obtain 

I Ve x (x),x - Vp{x) I = C(e 7 ) , 

for x G 5(0 , A — ce a ) with c > 0. We derive (2.5) by changing e\ by e in the previous 
estimate. Finally, writing 

p(x) _ (A + \x\) 

A dist(x, dV) A 

we get (2.8) for \x\ < A. 

D 

3 Rewriting the energy 

In this section, we prove equality (1.9), that is, the reformulation of the Gross-Pitaevskii 
energy of a two component condensate in (1.4), as the weighted Cahn-Hilliard energy for 
the pair (v, ip) defined by (1.7), plus the energy of the ground state r\ £ of a one component 
condensate. We start by giving the properties of the minimizers of S £ and the properties 
of the corresponding pairs (v £ , p £ ) defined by (1.7). 

Proposition 3.1. (i) Let {(wi, e , W2,e)}e>o be a sequence of minimizing pairs of £ £ in 7i 
satisfying (1.16). Then, each component is a non vanishing smooth function, and there 
is C > such that 

\\ui,s\\l°°(R 2 ) j ||w2,e||L°°(R2) < C (3.1) 

for every e > 0. Moreover, the pairs (v £ ,p £ ) are well defined by (1.7), verify the mass 
constraints (1.8) and we have 

(v s , p £ ) G Li Ploc (R 2 ■ (0, +oo) x [0, it]) (3.2) 

and 

sup ||i> e ||L°°(.K") < Ck for every K CC V . (3.3) 

£>0 

(ii) Conversely, let (v, tp) G Lip(M? ; (0, +oo) x [0, n]) satisfying (1.8) such that v , Vi> , V<p G 
L°°(M 2 ). Then, defining 

ui = i] £ v cos (y?/2) and -u 2 = Ve v sin ((p/2) , (3.4) 

we have (-Ui,-u 2 ) G H. and \ui\ 2 + \u2\ 2 > 0. 

Proof: (i) Let (ui, e ,ti2,e) be a minimizer of S £ in "H. Since £ £ (\ui t£ \, |ti2,e|) < £e(ui,e,U2,e), 
the pair of the absolute values satisfies the system 

- Au h£ + (e~ 2 V + £~ 2 u\ £ + g e u\ £ ) u h£ = X 1}£ u h£ (3.5) 

-Aw 2 , £ + [e~ 2 V + s~ 2 u\ £ + g £ u\ £ ) u 2 , £ = A 2 , £ u 2 , £ , (3.6) 
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where Ai j£ and A2, e are the Lagrange multipliers associated with (1.5). The strong maxi- 
mum principle yields that |ui j£ | and |«2, £ | are positive functions. Using standard elliptic 
regularity, we deduce further that Wi j£ and u 2)£ are non vanishing smooth functions. We 
use an argument in [20] to prove that Wi ;£ and u 2>£ are uniformly bounded in R 2 . Let us 
define w = £ _1 |wi i£ | — A £ . We have w G Lf oc (R 2 ) and Aw G L^IR 2 ). Kato's inequality 
and equation (3.5) give 

A(w + ) > sgn+(w) Aw > e~ 3 sgn + (w) ew (ew + eX 1 / 2 ) (ew + 2EX 1 / 2 ) > (w + ) 3 . 

Hence, — A(w + ) + (w + ) 3 < weakly in IR 2 and Lemma 2 in [12] yield w + < 0. We 
obtain |wi ;£ | < eX £ . Multiplying equation (3.5) by Ui j£ and then integrating we find 
A £ < 2 y/S e (\ui t£ \, \u2 t e\)- Since ( - Ui i£ ,'U2,e) verifies (1.16), from estimate (2.2), we derive 



< K, £ | < eJ (S £ (u 1>£ ,u 2 , £ ) - E £ (r] £ ^j + E £ (rj £ ) <C. 

We similarly prove that < |w 2 , £ | < C, so (3.1) is proved. Since r\ £ > 0, u\ >e and w 2)£ do 
not vanish in IR 2 , the pairs (v e ,(p e ) are well defined by (1.7) and v £ > 0. Since w l£ and 
M2,e are smooth, v £ and <p £ are locally Lipschitz functions so (3.2) holds. The definition of 
v £ and (1.5) give 

?7 2 t; 2 = «i + a 2 ■ 
From the definition of tp E , we infer that 



\u 1)£ \ 2 - \u 2 , £ \ 2 , , 

COS ^ = I ^^ 12 ' 3 - 7 



which, together with (1.5), yields 

T] 2 V 2 COS <f e = QL\ — a 2 ■ 

2 

Hence, (v £ ,ip £ ) satisfies (1.8). Finally, the estimate (2.4) gives r\ £ > ck > in K CC V, 
so (3.1) yields (3.3). 

(ii) Consider (v,ip) as in the statement and define (ui,u 2 ) by (3.4). Since (v,ip) verifies 
(1.8), relation (3.4) gives 



and 



\ui\ 2 + \u 2 \ 2 = / i] 2 v 2 = a 1 + a 2 



|2 I |2 111 1 

\U\\ — \u 2 \ = I r] £ v cos ip = «i — a 2 . 



Thus, (ui,u 2 ) verifies (1.5). We have |«i| 2 + \u 2 \ 2 > 0. Indeed, if it was not the case, since 
v > then tp should take simultaneously the values and it. Since v G L°°(]R 2 ), bounds 
(2.5) and (2.6) on r) £ give 



/ V (\ Ul \ 2 + \u 2 \ 2 ) <C f Vr] 2 <+oo. 

JR2 ' J R 2 
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We compute 

|Vtii| 2 <C7(t; 2 |V%| a + i7 e 2 |V%| 2 + A. 2 |V^ 2 ) • 

The right hand side of the inequality is integrable in R 2 because v , Vi> , Wip G L°°(IR 2 ) 
and r] £ G H 1 (R 2 ) C\ L°° (R 2 ) . Thus, u h£ G H^R 2 ). We prove similarly that u 2 , £ G H 1 ^ 2 ). 
We have proved that (ui t£ ,u 2t£ ) G "H. 

n 

We now prove the rewriting of the energy. 

Proposition 3.2. Let (ui,u 2 ) G "H satisfying \ui\ 2 + |^2 1 2 > 0. Defining (v,(p) by (1.7) 
we have 

€ e (u u u 2 ) = E £ (r] £ ) + F £ (v) + G £ (v, y?) , 

where E £ , F £ and G £ are given respectively by (1.1), (1.10) and (1.11). 

Proof: since |wi| 2 + |w 2 1 2 > 0, the pair (v,(p) is well defined. The definitions of v and (p 
yield 

\ui\ — rj £ v cos (ip/ 2) and \u 2 \ — r] £ v sin (ip/ 2) , (3.8) 

which give 



|2 , 1 |2 
Ul\ + \U 2 \ 


2 2 

= ^ 


1 |2| |2 
Pll l U 2| 


4 4 n 2 1 
= 4^ I 1 - COS V?} 


|4 , 1 |4 
Mil + |U2| 


= 2^ 4 ^ 1 + cos2v9 ^- 



(3.9) 



Since U\ and m 2 are real and do not change sign, we have |Vwi| 2 = | V|iii| | 2 and IVW2I 2 — 
|V|m2|| 2 - The relations in (3.8) give then 

\V Ul \ 2 + |V« 2 | 2 = |VH £ )| 2 + ^(v V£ ) 2 \V V \ 2 . (3.10) 

Replacing (3.9) and (3.10) in £ £ (ui,u 2 ) we get 



Se{u^U 2 ) = I j MVT] £ )\ 2 + ^Vny (3.11) 

+ \J \v 2 V 2 \V V \ 2 + ^ £ V{1 + COS^} + \geTiy{l ~ COS 2 ^} . 

The previous formulation of the energy is the one given by the spin formulation (see 
the introduction and [24]). We now show how the phase transition model is obtained. 
Performing an integration by parts, using (2.1) and the first mass constraint in (1.8), we 
obtain 



13 



\V{vr] e )\ 2 + ^V V 2 £ v 2 = v 2 r ]£ ^-Ar ]£ + ^Vr ]£ ^+r ] 2 \Vv\ 2 

= / v\ ( - Ar] £ + —Vr] £ + —7] 3 £ J - —rf £ v 2 + r] 2 £ \Vv 



2 



Ar Ion I Oi — iO -L 



f 2 JvWs+ JvlWv^-^y (3.12) 

Using again (2.1), together with the mass constraint for r/ e , we have that 

Ae 2(E £ ( V£ ) + ±- 2 J V t) . (3.13) 



S 2 



Replacing (3.13) in (3.12), and then (3.12) in (3.11) we get 
S £ ( Ul ,u 2 ) = E £ {r l£ )+ l -Jr 1 2 £ \Vv\ 2 + ^ 2 4{l-2v 2 } 



+ \ f \vWe\V<p\* + ^y{l + COSV} + \geVy{l ~ COSV} 



Completing the square for {1 — 1> 2 } we get 



£(ui )U2 ) = EM + ^frf^vf + ^il-v'Y-^y 



+ \ J \vW e \V<p\ 2 + ^^V{1 + cos 2 y?} + ^ £ 77 £ V{1 - cos 2 y?} 
Es{v) + \jvlWv\ 2 + ^ £ {l-v 2 } 2 



n 



+ ^ / ^\ 2 |Vy?| 2 + ^ £ 4 A £ (l - ^) {1 - cos 2 v?} , 
which finishes the proof. 

4 Upper bound inequality 

In this section, we consider the formulation of the problem in (v , </?) and call 

.F e (<;,<p)=F e (<;) + G e (t;,<p). 

We prove here the upper bound inequality for eT £ : 

Proposition 4.1. (Upper bound inequality for eJ^s) Let (p = 7t1a G X. There is 
a sequence of pairs (v £ ,ip e ) G Lip(M. 2 ; (0, 1] x [0, 7r]) ; converging as e — > to (l,y?) in 
L)ociV) x Lj^P), such that 

]i-msuveJF £ (v £ ,(p e ) < F(tp) . 
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The proof is based on Bouchite's paper [10], where he proves the r-convergence of an 
anisotropic phase transition Cahn-Hilliard energy. We point out that our weight r\ e de- 
pends on e and vanishes asymptotically on the boundary of T>. 

In a first step, we assume that (p = ttIa, where A is an open bounded subset of R 2 with 
smooth boundary such that 1-L l {dA fl dT>) = 0. Then, for any </? G X we approximate 
{ip = 7r} by this kind of sets. We conclude then thanks to a density argument. We remark 
that we do not consider here the mass constraints in (1.8). 

Before proving the upper bound, we recall some results about sets with smooth boundary, 
that can be found in Lemmas 3 and 4 of [26]. For an open set 4cK 2 with smooth, non 
empty compact boundary, let d be the signed distance to dA, defined by 

,, s _ J dist(x, dA) if x E A 
d \ x ) - | -dist(x,dA) if ie R 2 \A. 

For small t > 0, consider the neighborhood of dA given by 

N t = {xeM 2 ; \d{x)\ <t}, 
with boundary 

S t = {xeR 2 ; \d(x)\ =£}. 
For t > small enough, there is a diffeomorphism $ between N t and cL4x]0, t[ such that 

36 >0, det|£>$| > b. (4.1) 

We denote by $ the component of $ in dA. Moreover, d is a Lipschitz continuous function 
in N t and we have that 

\Dd\ = 1 a.e. in N t . (4.2) 

For small t > 0, define the measure 

lit^'LfDn s t ) . 

Notice that /i = % x L (£> n 9^4). As in Lemma 4 in [26], K 1 ^ n dV) = yields 



for every open Oct 2 , and 



liminf /i t (r2) > ^q{VL) 



Um^(P)=^ p). (4-3) 



Hence, as £ — >• 0, ji t converges weakly* to /io, which implies 



limsup / udfjit < / -urf/xo (4.4) 

for every upper semicontinuous function u : T> — > R with compact support (see Proposi- 
tions 1.62 and 1.80 in [5]). 

Denote ?? = y/p and for e > define f £ : M 2 x R x 1R 2 ->• R + by 
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fe(x,t, P ) = \V 2 AX)\P\ 2 + \vt(x){l -t 2 V- 

For \p\ = 1 and sGRwe also write f £ (x,t, s) = f £ (x,t, sp). 

The last step in the proof of Proposition 4.1 uses the following Lemma, which proof is 
given in the appendix. 

Lemma 4.2. Let A be a subset ofD with 1 A G BVi 0C (T>). There exists a sequence {Ak}km 
of open bounded subsets of M. 2 with smooth boundaries such that: 

(i) lim^oo C 2 ((A k n V)AA) = 0, 

(ii) limsupfc^, j v p 3 ' 2 \Dl Ak \ < f v p 3 / 2 \D1 A \, 

(ii) J A nV f i = J r A p and 'H}{dV fl dA^) = for k large enough. 

Proof of Proposition 4.1: we first assume that A is an open subset of 1R 2 with smooth, 
non empty compact boundary such that 

n 1 (dAndV) = 0. (4.5) 

(Step 1: construction of the pairs of test functions.) For T > 1, consider the approxima- 
tion of the optimal profile 

{tanh in (0, T) 
h in (T,T + i/t) 
1 in (T + i/t, +oo) 

where h is the unique cubic polynomial such that h(T) = tanh(T), h'(T) = tanh'(T), 
h(T + 1 /t) = 1 and h'(T + 1 /t) = 0. Computing explicitly the coefficients of h, we 
find that Wt is a nondecreasing function in R + , with uniform C 1 -bounds with respect 
to T G (1, oo). We extend wt to the whole real line by setting writ) = WT(—t) in R_. 
For e > 0, consider < m £ <^ e (to be chosen later), t £ = tanh - (m £ ) and define a 
modification of Wt near zero by 

_ / m £ in (0,t £ ) 
\ w T in (t e , oo) . 

Notice that w £ ^t has uniform Lipchitz bounds with respect to T G (1, oo) and e G [0, 1). 
We recall that T> = B(0, A) and we denote T>s = -8(0, A — 5). For y in dA\D$, we define 



and we write w T = Wq T . For small S > 0, we define R = Rs by 

r = {t+1/t) \Ij\- (46) 

Since % has uniform C 1 -bounds with respect to T G (1, oo), while p is a smooth function 
in T>s, for every y G dA D 2^,5, there is an open neighborhood S of y in cM fl P,5 such that 
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R rR 

fo(x, w y T (t), (w y T )'(t)) dt< f (x, w^(t), («#)'(*)) dt + 5 ViGE, VT > 1 . 

Jo 

Hence, thanks to the compactness of dA D T)$, there is a finite family {Sj}^ 1 of open 
disjoint subsets of dA fl V$, and a corresponding family of points y« G Ej, such that 

H 1 (dAnV s \\JVi\ =0 (4.7) 

and 

/ f {x,w*(t),(w*)'(t))dt< [ fo(x,w* T (t),«)'(t))dt + 5, (4.8) 

Jo Jo 

for every x G E„ T > 1 and 1 < % < N. We will use the functions wf T to define the first 
test function, so we have to interpolate between the different Sj's. Define first So = dA\T>s 
and y = (\-5, 0) G dV 5 . For small £ > define S^jiGS,; dist(x, <9E;) > £}. Clearly, 

^(E^Si) -^ as £^0. 
In particular, we can take £ = is such that 

RH 1 piXEt) = o s ^o(l) (4.9) 

for every < i < N. Consider then {^}il such that & e C°°(9A, [0, 1]), 

N 

^0, = 1 on dA and £ = 1 in Ef . (4.10) 

We deduce a smooth partition of the unity on N t by setting 6i = 6i o <E> and we define 



1 in M 2 \JV e 



etf 



Z£oft(x)< r (^) in iV efl 



Since w^ T is a nondecreasing function, while p is a radial decreasing function, (2.8) and 
the fact that dist(?/j,X>) > <5 yield 



<tU. b = <AR) = Ve,T l(T + 1/T) ^^ j > «, e>r (T + 1/T) = 1 , (4.11) 

so f £ is a continuous function. Moreover, since w £j t has uniform Lipschitz bounds with 
respect to T G (1, oo) and e G (0, 1), there is C > such that for e small enough, 

C 
IKIIc°.i(iR 2 ) < — • (4.12) 

We also define 
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71 if x e A\N ei 

Ve(x) = { C( d WM) if xeN £ie , (4.13) 

if x E R 2 \(A U N ei J 

where £(t) = */z(l + *) and t £ = ( 2 /\) 1/2 t £ . We clearly have that (y e ,<p e ) E Lip(R 2 ; (0, 1] x 
[0, 7r]), and that (v £ ,(p e ) converges as e — > to (l,<p) in L\ oc (V) x L\ oc (V). 

(Step 2: estimating the energy eG £ .) The function (p £ is constant in R 2 \N £ i , so G £ (v £ , y? £ ) = 
G £ (v £ , ip £ ; N e i s ). Since w £> t is a nondecreasing function while p has a global maximum at 
zero, for every x G N E f we have 

( A |d(s)|\ , , 

= u> £ , T ( -^ i -j ii J < w e , T (t e ) = m £ , 



<r( W f)<<ri 


^|d(x)| 


so t> £ < m £ in iV £ j £ . Hence, 





G £ (v £ ,tp £ ) 



<o / ^e|V^ e | 2 + ^m^ e {l-cos 2 (v? £ )}. 



Then, the definitions of <p £ and i? £ , together with the fact that i] £ is uniformly bounded, 
yield 

G £ (v £ , Ve ) < C {m 2 £ {et £ )- 2 + g £ mf) C 2 {N de ) . 

Using (4.1), we have 



{et £ )-'C\N de ) < C(et £ )- L / dt |det (£>$) _1 | dU 

J-ei E JdA 

< b- l H\3A) 

= eWo(i)- 

For e small, we have i £ = ( 2 / A )^ 2 tanh~ 1 (m £ ) = O £ ^o(m £ ). Hence, 

G £ (v £ , ip £ ) < C {m^' 1 + m\ g £ e) . 
Taking m £ = (g £ e 2 )~^ 4 , G £ (v £ ,(p £ ) < Cg £ e~ 3 ^ 2 , and after (1.6) we obtain 

eG e (v e ,<p e ) < C(g £ s 2 )-^ = o £ ^ (l) . (AAA) 

(Step 3: computing the energy eF £ .) Since v £ is constant out of N £ r, we have 



eF £ (v £ )= / £ + / £ + / £ . (4.15) 

iN ER nv 5 JN eR n(v\v s ) JN eR \v 



where 



4> £ {x) = -f £ (x,v £ ,eDv £ ). 

Considering (4.12), for e small enough there is C > such that |Vt> £ | < c /e. Estimates 
(2.4), (2.7) and (2.8) thus yield 
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/ £ < c(max{p + p 2 } + C s e 2 \lne\)e- 1 C 2 (N £R n(V\V s )) 

JN sR r\{v\v 5 ) v dv t J 

< C(8 + e 2 \\ne\)e- 1 C 2 {N eR ). 



Using (4.1) we have 



e- l C 2 {N eR ) < Ce~ l [ dt f |det (D®)' 1 ] dU l 

J-eR JdA 

< CRb- l U\dA), 
so (4.6) yields 

lim f <P £ <C5Rb- 1 y}{dA) = o s ^ Q {l). (4.16) 



'N eR n(v\v s ) 
Similarly, using (2.6) we have 



/ <Pe < C svv(r 1 2 £ +4)Rb- l U l {dA) 

Jn eR \v m?\v 



< Ce^Rb^H^dA) (4.17) 

= o e _,o(l) • 



Now, remember the interpolation from (4.10). We have 

JN £R nv s ■ , \JN eR nBi JN eR nc t / 



'N eR nv s - =1 \JN eR nBi JN eR nc t 

where 

^ = {x 6 P tf ; $(s) G Ef} and Q = {i6D f ; $(s) G E;\Ef} . 

As before, since 7/ e is uniformly bounded in M 2 with respect to e G (0, 1), (4.1) yields 

/ <f>e < C- [ dt I \det(D^>)- l \dV} 

JN sR nCi £ J-eR ./£i\Ef 

< cRb^n 1 ^^). 

Hence, after (4.6) and (4.9), we obtain 

£ = o^ o (l) (4.1? 



'N eR nd 
for every e G (0, 1). In N sR H -B; we have v £ (x) = wf T (\ d ( x )\/e). Using (4.2) we write 

£ = I \D{\d\/e)\{x)f £ {x,wf T {W*)\/s)i « T )'(W*)IA)) • 

AT ei jnBi JN eR nBi 
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The coarea formula from Proposition 1.5 yields 

-R 



dt / f e {x,wf T {t), (wf T y(t)) dfx £t (x) . 

N eR <lBi J-R JT> 6 nB, 

We thus have, 

e< [ dt f f (x,w^(t),(w^)'(t))d^ t (x) + Ri + Ri. (4.19) 

N eR nBi J-R Jv s nBi 

The first error here before comes from the modification of wt near 0. Using (4.3) and the 
definition of t £ we compute 



R\ < CAl——t £ sup \\nt\\ 
PKVi) te(o,sR) 

- c \ j\ t£ su p ll/ 1 *!! 

V OA te(0,eR) 

= o e _ (l). 

The second error appears when replacing f e by /o, so using estimates (2.4) and (2.8) 
together with y { E T>$, there is C$ > such that 



R 2 e < C s e 2 1 hie | R sup t6(0)efl) ||/i t || = o £ _> (l) • 
Using Fubini's formula, we rewrite (4.19) as 



•e ^ / J-CsOB 



< / l-DrnR,fx) 



/ (x,<(t), «)'(t)) dt j dHet(x) + O e _ +0 (l) . 



Af £fl nB, Jv \ J-R J 

The set T>$ fl Sj is close and the inner integral is a continuous function of x. Hence, the 
function inside the outer integral is upper semicontinuous function of x. Inequality (4.4) 
thus yields 

limsup/ (j) £ < I \lv s nBi(x) f (x,w%(t),(w%)'(t))dt) d/J, (x) 

s->0 jN eR C\Bi Jv \ J-R J 

R x 

f (x,w$(t), (w^)'(t)) dt dno(x) . 
v s r\T,\ \J-R J 

Notice that since yU is supported in dA, we replaced Bi by Vs fl Sf . From (4.8) and since 
Wt is an even function, we have 

limsup / e < 2 / I / f (x,w T (t),(w T )'(t))dt + 8 ) d(j, Q (x) 

e^O JN eR nBi Jv s C\Y>l \Jo J 

< 2 [ ( [ f (x,w T (t), (w T )'(t)) dt + 6] dn (x) . (4.20) 
Jvscfcl \Jo J 
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(Step Jf.: upper bound) Putting together (4.14)-(4.18) and (4.20) we get 



limsupej;^,^) < 2 V / ( 



fo(x, u%{t), (i4)'(t)) dt + S] dno(x) + o^o(l) 



Now, we take a sequence T = T$ such that T^ — >■ oo as 5 — > (notice that (4.11) still 
holds). Then, Fubini's formula and dominated convergence theorem, together with (4.7) 
and (4.9), yield 

lim (limsupe.F e (t; e ,^ e )J < I f (x, w x (t), (w x )'(t)) dt j dfi (x) . 

Remembering the definitions of /o, fJ,o an d w 1 '■> (1-14) yields 

lim I limsupeJ-' e (f e , v? £ ) J < 2a l p 3 ' 2 (x)d / H 1 (x) . 

We conclude thanks to a diagonal argument (see Corollary 1.16 in [7]): there exists a 
sequence S £ — ?- £ ^o 0, such that as e — > 0, (v e ,6 e ,<Pe,5 e ) converges in L\ oc {T>) x L\ oc {V) to 
(1, ip), and 

limsupeF^SeivejctPeje) < 2o " / p 3/2 (x)dn 1 (x). 

e^O JVndA 

(Step 5: approximation of A by Cacciopoli sets) We end the proof using Lemma 4.2, the 
proof of which is given in the appendix. We remove the condition (4.5) and we only 
assume that A is a set with locally finite perimeter in T>. Consider <p k = fflA k with 
{^4fc}fceN the sequence from Lemma 4.2. From (i), (l,y2 fc ) converges to (l,ty?) in LMj)) 
and we have j A p = a 2 f° r k large enough. Hence, from steps (l)-(4), there is a sequence 
(v k ,tp k ) -» (lj k ) in L\ oc {V) x L\ oc {V) such that 

\uneF e (v k e ,ip k ) = F{ip k ) . 
Using (ii) from Lemma 4.2 we obtain 

lim sup \Y\m.eJ : e{vl,ip k ) ) < .F(</?) . 

As in step (4), we conclude thanks to a diagonal argument. 

□ 



5 Lower bound inequality and compactness 

The proofs in this section are based on geometric measure theory techniques. We make 
the lower bound on lines and then use the slicing method, which can be found in [6] or 
[11]. The last part of the proof of the lower bound is inspired by the ideas in [4]. 
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5.1 Lower bound on lines 

Consider an open set A CC T> and let v G S 1 be a fixed direction. We call tc u the hyper- 
plane orthogonal to v, and A v the projection of A on tt v . We define the one dimensional 
slices of A, indexed by x G A v , as 

A x = {t e R; x + tv e A} . 

For every function / in D, we define f x as the restriction of / to the slice A x , defined by 
fx(t) — f(x + tu). For (v,<p) : A x — > (0, 1] x (0, re), we define the energies 

Fe(v,A x ) = l[vUVv\ 2 + ±- 2V l x {l-v*}\ 

G £ (v,<p;A x ) = - r] 2 x v 2 \Vip\ 2 + T]l x v A g £ {l-cos 2 (ip)} and 
o Ja x 

T £ (v,<p;A x ) = F £ (v;A x ) + G £ (v,ip;A x ). 
Similarly, for p G BV(A X ; {0, 7r}) we define 

HV,A X ) = - f pfd\D<p\. 



A, 



With the previous notations, we have the following result: 
Proposition 5.1. Let (v £ ,tp e ) G Lip(A x ; (0, +oo) x [0, 7r]) such that 



sup||w £ || L oo (Aa;) < oo (5.1) 

£>0 



and 



eF £ (v £ , (p £ ; A x ) < oo . (5.2) 

Then, there is tp G SBV(A X ; {0, 7r}) such that 

(v £ , <p £ ) -> (1, </?) in L\A X ) x L\A X ) (5.3) 

and 

limini eJ r £ (v £ ,p £ ; A x ) > T(p ; A x ) . (5.4) 

Proof: (Step 1) Using that A CCV and estimate (2.4), there are C\, c 2 > such that 

i] 2 x > p x - cie 2 \lne\> c 2 in A x . (5.5) 

Hence, the definition of T £ {; ; A x ) and (5.2) give 

/ |1 -v £ \ < — ^-e 2 T £ {v £ ,p £ \A x ) = o e _K)(l), 
so f e — )■ 1 in L 1 {A X ). Similarly, after (5.1) v £ < C in Aj, so (1.6) yields 
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/ V £ \l- COS 2 (ip £ )\ < 2-Fe(v e ,<p e ;A x ) = O e _M)(l). 

JA X 9e C 2 

Hence, up to a (not relabeled) subsequence, ip £ —¥ <p a.e. in A x , with <p : A x — > {0,n}. 
This, together with A x CCV, gives ip £ — > ip in L 1 (A ;C ). We have proved (5.3). 

(Step 2) We now prove the lower bound for the energy. Let to G Sip. For 5 > define 

J 5 = A x n (t - 5, t + 5) , 
and suppose that 

inf (inf^(t)} >c 3 >0. (5.6) 

teJa le>0 J 

Then, for every e > and every t G Jg, v e (t) > C3. Hence, using (5.2), (5.5) and the 
coarea formula (1.21), there is C > such that 



(c 2 l cf + cfc^)— 2 > [ |Vy.| 2 + {l-cos 2 (y.)}> [* dt {l-cof®} 1 * [ \DlwJ, 

9e£ Jj s Jo Jj s 

(5.7) 
where W £>t = {t G A x ; (p s (x) < t}. Since <p £ converges to ip a.e. in A x , we get 

% £ , ( -> 1{<^=o} in ^(AJ, 

for a.e. t G (0, 7r). Hence, the lower semicontinuity of the BV norm with respect to the 
/^-convergence, together with (1.6), (5.7) and Fatou's lemma, gives 

> / dt{\ -cos 2 (t)} 1/2 liminf / \Dl Wet \ 
Jo £ ^° Jjs 



'Js 
L {v=o}|- 

* -JJs 

Thus, 



7T / 



o = ^ (J,n 1 s V 9)>^ ({t }) = i. 

This contradiction implies that (5.6) can not be satisfied. We derive that for every 8 > 0, 
we may extract a subsequence (not relabeled), such that exists {t £ } £>0 C J 5 with 

t £ ->■ t G J«5 and v e (i e ) ->■ as £ ->■ . (5.8) 

(Step 5; For e > 0, define if = {t G J s ; ±(t £ - t) < 0} and v± : J 5 -> (0, 1] by 

The definition of J- £ , estimate (5.5) and the fact that vf is constant in If while equal to 
v £ in If yield 



V2sF £ (v £ , V£ ;J s ) > J p?(\(v+)'\ |1 - «) 



+ ^ 2 ' M(^-yiu-(^-) 2 in-o^n(i). 



e y 11- ve / 1 ; ~e-»-0 
Js 
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Using the coarea formula (1.21) we obtain 

V2eF e (v e ,<p e ;J s )> I dt(l-t 2 ) [ pf (\D1 V+ | + \D1 V - |) + o £ ^ (l) , 

JO Jjg V e ' e ' J 

where V^ t = {t G J$; vf < t}. Since t £ — ¥ i , v £ (t £ ) — > and v £ (t) — >■ 1 a.e. in 
J<5, l y ± — > 1/=f in L x (Js), where J 1 * 1 = {t G Ja ; ±(t — t) < 0}. Hence, the lower 
semicontinuity of the BV norm with respect to the /^-convergence and Fatou's lemma 
give 

]imini£T £ (v £ ,ip £ ;J s )>^= f dt{l-t 2 ) j pf (|Dlj-| + plj+l) = a2pf (f ) ■ (5.9) 

Moreover, since p^ > C4 > in Ar, we have 

liminf sT £ (v £ , (p £ ; J5) > 2crc 4 > . (5.10) 

(Step 4) Let T = {£1, . . . , t n }, n G N, be any finite subset of Sip. For i G {0, 1, . . . , n} we 
define 

J\ = A x n (*i - <J, ii + «J) . 
Consider <5' > such that JJ fl J^ = for i 7^ j and let S G (0, <5'). From (5.10), we have 

— liminf £.F £ (?; £ ,<p £ ; A x ) > n. 

2a e^o 

Therefore, using (5.2) we derive that n is bounded, so Sip is a finite set and p G SBV(A X ). 

(Step 5) Finally, write Sip = {t±, • • • ,t N }, N G N. Reasoning as before, for 5' small 
enough and S G (0,6'), (5.9) gives 



liminf eT £ (v £ ,p £ ;A x ) > 2cr V ]pf(U) , with U G J\ 



j=l 



Since p x is a continuous function, taking the limit 5 — > in the previous inequality we 
obtain 

N 

\immie7 £ {v £ ,p £ ;A x )>2ay^pf{t i )=2a / pf d\D<p\ . 
We have proved (5.4). 

n 



5.2 The slicing method 

Using the slicing method, we now prove the compactness and the lower bound inequality 
for eJ- £ . 
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Proposition 5.2. (Lower bound inequality and compactness for sJ- £ ) Let (v £ , tp £ ) G 
Lipi oc (R 2 ; (0,+oo) x [0,7r]) snc/i £/ia£ 

sup 1 1 f e || £,«>(#) < Ck for K CC V (5-H) 

£>0 

and 

mpsF £ (v £ ,tp £ ) < oo. (5.12) 

£>0 

Then, there is tp G X such that 

(v £ , tp £ ) -> (1, </?) m LLp) x LLP) (5.13) 

and 

liminf sT £ {v £ ,tp £ ) > T(tp) . (5.14) 

Proof: arguing as in the proof of (5.3) of Proposition 5.1, there exists tp : T> — > {0, 7r} 
such that (5.13) is satisfied. Consider an open set A CC T>, and fix u G S l . For x G A„, 
we define (v £yX , tp £jX ) : A x ->■ (0, 1] x (0, 7r) by 

{Ve,x, Ve,x)(t) = (v £ , tp £ )(x + tu) . 

Since A CC T> and since i> £ , r^ £ are non vanishing continuous functions, for fixed e > 
(5.12) yields 



C>£j;(^,^;A) > £ -mi V 2 j \Vv £ \ 2 + £ -miv 2 V 2 j \Vtp £ \ 2 



> c £ , A U \Vv £ \ 2 + J |Vy? £ | 2 r . 
so v £ and tp £ belong to W 1,2 (A). Hence (see [15], Section 4.9.2), 

v' £X (t) = D u v £ (x + tu) and f' e ,x{t) = D v <p e (x + tu) 

for a.e. t G A x , for £ x - a.e. x G A„. Using then |Vt> £ | 2 > \D u v £ \ 2 , we get the slicing 
inequality 



eJ r £ {v £ , tp £ ■ A) > / eJ r £ {v £tX , tp £>x ; AJ cb . 



(5.15) 



From (5.12), for £ 1 -a.e. x G A u , sJ- £ (v £}X ,tp £tX ;A x ) is uniformly bounded with respect to 
e. Thus, after Proposition 5.1, for £ 1 -a.e. x G A v there is tp x G BV(A X ; {0, 7r}) such that 

(v £ , x , <p £ , x ) -> (1, y x ) in L^xi 1 ^) (5.16) 

and 

liminfeJ r e (u eil! ,y> ea .;A I .) > J 7 ^ ; A„) . (5.17) 

£>0 
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The function p denned in (5.13) is the L l (A) limit of p E , so for £ 1 -a.e. x G A u , tp x 
coincide with the restriction of tp to A x . Therefore, since the vector v is taken arbitrarily, 
tp G BV(A) (see Proposition 6.9 in [4]), and since A is any open relatively compact subset 
of T>, we derive that tp G BV\ oc {T>\ Using (5.15), (5.17) Fatou's lemma and Fubini's 
formula, we also obtain 



\m\'vnie7 e {v e ,p e ) > / T(l,tp x ;A x )dx 

£ ^° Ja v 

2a 



K J A 

2a 



dx / p z J. 2 d\Dtp\ 
Ja x 

[ pfdiC^A^lDtpxlLAx). (5.18) 

Jv 



Now, for every e > 0, let p £ be the energy distribution in T> associated with the pair 
(v £ , tp £ ), that is, the positive Radon measure which for every Borel set E C M. 2 is given by 

f x e (E)=eT e (v et (p e ;EnV). 

From (5.12), the total mass ||/x e || is uniformly bounded. Dc La Vallee Poussin compactness 
criterion (see [5] , page 26) gives then that (up to a subsequence) p £ converges weakly* to 
some finite measure \x on T>. We claim that 

fi > 2ap i/2 ■U 1 \-Sip. 

We will prove this using Besicovitch derivation Theorem (see [5], page 54). First, after 
(5.12) for every K CC V there is R K G (0, A) such that 

< fi(K) < fi(B(0,R K )) < ]imin£ ii e (B(0,R K )) < oo. (5.19) 

Hence, \x is a positive Radon measure in T>, and for 7i l -a.e. x G Sip the limit 

f(x) = lim —j&^LL - (5.20) 

Jy ' r^Q+'H 1 {B r (x)r\Sp) ^ ' 

exists, and we have 

pyf-U^Sip. (5.21) 



Let Xq G Sip fl A. Since A GG V, B t {xq) G T> for r small enough. We assume 1 moreover 
that fj,(dB(x ,r)) = 0. Proposition 1.62 in [5] and estimate (5.18) yield 



n(B r (x )) = lim n E (B r (x )) 



> — / pfd(C l ^A u ®\Dp x \^A x ) 

71 JBJxo) 



(X ) 

> ^ inf pf [ d{C 1 \^A v ^\D(p x \[.A x ). (5.22) 

"V(xo) is r (x) 



7T B r ^xo; JB r (x) 

^n fact this holds for all r except countably many (see [5], page 29). 
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In Proposition 5.1 we proved that for £ 1 -a.e. x G A u , tp x G SBV(A X ) fl L°°(A), and 

/ drr / n°(S Vx ) <oo. 

J A v J A x 

Hence, after Theorem 2.3 in [6], <p G SW(A) n L°°(A) and 

/ d(£ 1 LA I/ (8)|D ¥ > x |L^ x )=7r / K^,//)!^ 1 , (5.23) 

J B r {xo) J B r (xo) 

where v^ is the measure theoretic inner normal to the Caccioppoli set {tp = tt}. Putting 
(5.23) in (5.22) we obtain 



fi(B r (xo)) > 2a inf pf [ \(v^v)\dU l 

B r-(*o) JB r (x ) 



> 2(7 inf pf inf K.^JjtflB^n^). 

B r (x ) B r (a;o)n5</3 



Since SV is a rectifiable set in A, for T^-a.e. xo G S<p, v v (x) is continuous in B r (xo) for r 
small enough. Thus, taking v = v^ixo) and since p is continuous, we get 

p(B r (x )) 3/ 

r^0+ fi 1 {B r {Xo) n £</?) 

for T-L l -a,.e. xq G SV>- Hence, (5.20) and (5.21) yield the claim. The definition of /x gives 
then 

liminf e.F e (t; e , p £ ; A) = liminf /i e (A) > /i(A) > 2ct / p^d^ 1 . 
Finally, taking an increasing sequence {A^eN with Ak CC V, we get 



liminf ej;(t; e ,v9 e ) > 1o I p A,2 (m} , 

£ ^° JSwnv 



which gives (5.14). 



□ 



5.3 Remark about a lower bound for v £ in the transition zone 

We end this section with a discussion about the infinimum of v £ in the transition zone. 
Let {(v £ , <p £ )}e>o be a sequence of minimizers of T e , and let p G BV\ oc (V ; {0, n}) be the 
Lj^-limit of tp e given in (1.17). Let K CC T> be an open smooth set, with non negligible 
intersection with Sip, that is, 

u^KnSip) >o. 

For every e > 0, we define 

m £yK = inf v £ (x) . 

x£K 
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We would like to obtain an upper bound for m £ ,K, in connection with an open question 
in [8], namely 

me,K<C K (g £ e 2 )-^. (5.24) 

If we assume that we have the upper and lower inequalities for each s > 0, that is 

eF £ (v e ) < T(<p) (5.25) 



and 



eF £ {v £ ) > F(<p) , (5.26) 



we can give estimates on G £ in order to obtain the upper bound for m £ K- So assume that 
we have (5.25) and (5.26). On the one hand, estimates (2.4) and (2.8) give then 

eG £ (v £ ,ip £ ; K) > -y/g^m\ iK (inf p i/2 - C K e 2 \ \ne 2 \) / |Vy? e |siny9 £ 
> C K \/g £ £ 2 m% K / |Vy? e |siny9 e . 

JK 

We claim that the integral here below is bounded away from zero. Indeed, if this not the 
case, we will have 

liminf / I Vy? e I sin ip e = . 
£ ^° Jk 

Hence, since <p £ — > <p in ^(K), the coarea formula together with the lower semi continuity 
of the BV norm imply the contradiction 

pn r pn p 

= liminf / sintdt \Dl {v>c<t} \ > / sintdt \D1 W=0} \ = 2V}(S^ n K) . 
e ~ >0 Jo Jk Jo Jk 

We thus derive that there is C' K > such that 



eG £ (v £ , <p £ ;K)> C' K ^g^ m% K . (5.27) 

In the other hand, by inspection of the proof of Proposition 4.1 (see estimate (4.14)), we 
see that the pair of test function (v £ , ip £ ) satisfies 

eG £ (v £ ^ £ ;K)<C(g £ e 2 )-V\ (5.28) 

Hence, considering (5.25)-(5.28), together with the fact that (v £ ,ip £ ) minimizes J- £ , we 
obtain 

2a f p 3/2 \D<p\ + C K y/g^mlje < eT £ (v £ , Ve ; K) 
Jk 

< eT £ (v £ ,<p £ ; K) 

< 2a f p^\D^\+C(g £ e 2 )-y\ 

Jk 
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Multiplying both sides of the previous inequality by (g £ e 2 ) ' 4 we find the upper bound 
(5.24) for m\ K . 

However, we are not able to prove (5.25) and (5.26) as such because of the error terms. 
Indeed, the proof of the upper bound of Theorem 1.1 says that there is a sequence 
{(v £ ,p £ )} £>0 such that 

limsupeF £ ({} £ ) < F(ip) . (5.29) 

In the proof of (5.29), we first approximate the locally Cacciopoli set A — {ip — n} by 
characteristics functions of open sets Aj, with compact smooth boundary. This gives a 
small error in terms of k G N in the upper bound inequality (5.29). Then for each k G N, 
we construct a test function for which (5.29) holds, up to a small error term depending 
on a parameter 5 > 0. In these two steps, we use diagonal extraction arguments in order 
to get rid of the error terms, so it is not possible to compute them explicitly. Similarly, in 
the proof of the lower bound of Theorem 1.1, we use the compactness of bounded Radon 
measures, so we cannot estimate the error term in the lower bound inequality 

liminf eFJv £ ) > F((p) . (5.30) 

6 Proof of the T-convergence for e (S £ (-) — E e (rj e )) 

6.1 Proof of the compactness and the lower bound inequality in 
Theorem 1.1: 

let {(wi, e , u 2t£ )}e>o De a sequence of minimizers of S £ in H satisfying (1.16). From Proposi- 
tion 3.1(i), the pairs (v £ ,<p £ ) are well defined by (1.7), belong to Lipi oc {R 2 ; (0, +oo)x [0, n]) 
and satisfy (5.11). Proposition 3.2 yields eT £ {y £ , <p £ ) < oo. Thus, the hypotheses of 
Proposition 5.2 are fulfilled by (v £ , (p £ ) and we have 



{v e ,<p e )^(l,<p) in L} oc {V) x L} oc {V) 

with ip G X, and 

lim inf e T £ (v £ ,ip £ ) > T( p ) . 
Equality (1.9) yields then 

liminf £ (£ £ (u li£ ,u 2 , £ ) - E £ (r] £ )) > F{p) . 
Finally, using identity (3.4) we get 

(ui, £ , u 2 , £ ) -+ ^TP (!{v=o}, 1-W=n}) in L} oc (V) x L) oc (V) . 



D 



In order to prove the upper bound we have to work a little more. We first modify the 
pairs of test functions from Proposition 4.1 to make them satisfy the mass constraints 
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(1.8). We prove then that this modification do not change the limit of the energy. We 
finish by verifying that the pairs of modified test functions are the image by (1.7) of a 
pair in "H, and we conclude using Proposition 4.1. 



6.2 Proof of the upper bound inequality in Theorem 1.1: 

(Step 1 : Modification of the pairs of test functions) With the notations from the proof of 
Proposition 4.1, we write N £ = N £ r 6 and we define (v e ,(p e ) the sequence of pairs of test 
functions such that 

\imsupeJ 7 £ (v £ ,(p £ ) < F(<p) . (6.1) 

Consider k G C°°(1R + ; [0, 1]) with supp k C (0, 1) and k — 1 in (0, l /2). Since A is a non 
empty open set, there is B = B ro (x ) CC AnV. For £ G [— 1, 1] and r G (1/2, 1), define 
«e = K e,e,r by 

K £ (x) = E T £K,(\ x - x o\/r ) . 

We define then v £ = v £ + K s and v £ = c £ v £ , with c £ = ||^ £ '0 e ||^ 2 . For e small enough N £ 
and B are disjoints. We estimate 



c- 1 = 1 + / ^ e 2 - 1) 

= 1 + 2 f V 2 £ K £ + f r, 2 M - 1) + f v y £ 

J B J N E JB 

= 1 + 2 f rf £ K £ + 0{e) + 0(e 2r ) . 

JBo 

Hence, using that r G ( 1 /2, 1) we get c 2 = 1 — r £ with 

Bo 



r £ = A 77X + (9(e) = 0(e T ) . (6.2) 

•/Bo 

Notice that for e small enough, r £ may be positive or negative depending on the sign of L 



The definition oiw y eT insures that v £ > 0. The first mass constraint in (1.8) is immediately 
satisfied by the definition of v £ . Remember the definition of (p £ in (4.13). For the second 
mass constraint we write 



c £ 2 / rf £ vl cos cp £ = / Vs(1ri\(aun e ) ~ 1a\(n e ub ) + 1n e ub v £ cos <p e ) . 

JR2 JIR2 

Adding and removing 1jv s \a?7 2 , lAf e uA^ 2 and 1b V 2 m the previous integral, we get 

cf f r,y £ cos^ = / ?? 2 (1 R2 - 21 A ) + / V 2 (v £ + k £ ) 2 - 1 (6.3) 

^R 2 Jm. 2 Jb 

+ / i] 2 £ (v 2 £ cos Lp £ -1 A + 1r2\ A ) . (6.4) 
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For the third term in (6.3), we have that r] £ , v £ and cosy? e are bounded while C 2 (N £ 
(9(e). Hence, 



r] £ (v £ cos ife - 1 A + 1 R 2\ A ) = 0(e) . (6.5) 



N, 



For the first term in (6.3), using that J R2 rj 2 = 1 = a.\ + a 2 and that § VnA p = oc 2 , we 
obtain 



/ V 2 e {l R 2 - 21 A ) = e*i - a 2 + / (rf £ - p) + f 

JR 2 JAnV JA\D 

Using (2.5) we get, for a G (Y2, 3 /5) and 7 G (Y2, 3 /4), 



^ 2 - 



(%-p) = / (vi-p)+ / (^ 2 -p) 

An© ^AnB(0,A-e«) J (AnV)\B(0,X-e a ) 

= 0{e r ) + 0{e a ). (6.6) 

Moreover, from (2.7), we have 77^(0;) < i] 2 (x a ) in A\D, with x a G dB(0, A — e a ). From 
(2.5) and (2.8) we get 

if £ (x) < rjl(x a ) = rjl(x a ) - p(x a ) + p(x a ) = 0(e a ) , 
so using that A is a bounded set we obtain 

Vs = 0(e a ) . (6.7) 

A\V 

For the second term in (6.3), the definitions of K e and r £ yield 

77 £ 2 « £ (2 + k £ ) = l -r £ + 0(e) + 0(e 2r ) . (6.8) 

Putting (6.5)-(6.8) in (6.3) and considering (6.2) we get 

/ r] 2 v 2 cos ife = a>i - a 2 + -r £ + 0(e 13 ) , 

JR 2 ^ 



c; 2 



where (3 = min{l, a, 7, 2r} = min{a,7} G ( 1 /2, 3 /b). Hence, (6.2) gives 

7] 2 v 2 cos ip £ - (ai - a 2 ) = ( ~ _ ( a i _ "2) J r e + 0(e p ) . 

Suppose now, without loss of generality, that a.\ — a 2 < l J2. The definition of r £ and k £ , 
together with (2.4), (2.8) and B CC Af]V, give then 
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|r e | > 4inf7? 2 / k 2 £ + 0(e 



I 4 

D ° Jb„ 

> 4mi V 2 \£\e^ [ 4 + 0(e) 

> c\£\e T + 0(e), 

for some c > not depending on e. Hence, if we take £ = 1 in the definition of k £ , for e 
small enough we have 



r\ £ v £ cos (p £ — («i — a 2 ) > c e T (l + £ 
Analogously, taking now £ = — 1, we get 

/ rgv* cos^ e - («x - a 2 ) < cV(-l + e 1 "" + e p ~ T ) 
Since /3 G ( 1 /2, 3 /5), we can choose r G (Y2, /3) and obtain 

rj 2 v 2 cos <^ e > ol\ — a 2 if £ = 1 



and 



i] £ v £ cosip £ < ai — a 2 if £ = —1 . 



Hence, there exists £ £ G (— 1, 1) such that for e small enough, the associated pair (v £ , ip £ ) 
satisfy the second mass constraint in (1.8). 

(Step 2 : Computing the energy). We now compute the energy of (v £ , (p £ ). We recall that 
N £ i is the transition zone of <p £ defined in (4.13). For the energy G £ , we have that tp £ is 
constant out of N £ ^, while v £ = c £ v £ in N e f s with c £ = 1 + 0(e T ). Hence, 

eG e (v e ,<p e ) = (1 + 0(e T ))eG e (y e ,<p e ) . (6.9) 

For the energy F £ , we have that v £ = c e (l + K £ ) in B . The definition of k £ gives then, 
\Vv £ \ 2 = 0(e 2T ) and {1 - v 2 } 2 = 0(e 2r ). Hence, 

eF £ (v £ ; B ) = 0{e 2r - 1 ) = o e - (l) • (6.10) 

In M 2 \(iv; n S ), we have that v £ = c £ , so |Vu e | = and {1 - v 2 } 2 = 0(e 2t ). As before 
we get 

eF £ (v £ ; R 2 \(N £ n S )) = ©(e 27- " 1 ) = o e -o(l) • (6.11) 

In JV £ , we have that t> £ = c e -u e . Hence, |Vw e | 2 = (1 + 0(e T ))\Vv £ \ 2 and {1 — v 2 } 2 = 
(1 + 0(£ T )){1 - v 2 } 2 + (D(e T ), which gives 

eF £ (v £ ; JV e ) = (1 + 0(e T )) eF £ (v £ ; JV e ) + 0{e T )e^C 2 {N £ ) 

= (l + O(e r ))eF £ (v £ ;N £ ) + o £ ^ (l). (6.12) 
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Since v £ is constant out of N £ , we have F £ {v £ ) = F £ (v £ ; N £ ). Putting together (6.1) and 
(6.9)-(6.12), we obtain 

limsupeJ^Ue, ip £ ) = UmsupeJ : ' £ (v £ , ip £ ) < F(tp) . (6.13) 

(Step 3 : identification of (v £ , <p £ )) The pairs of test functions satisfies the hypothesis 
from Proposition 3.1(ii), so defining (tii j£ ,ti2,e) by (3.4) we have (ui i£ ,u 2j£ ) £ H and 
u ie + M 2e > 0- Hence, after Proposition (3.2) relation (1.9) holds, and (6.13) yield 



limsupe (£ £ («i, £ , w 2 , £ ) - E £ (rj £ )) = \imsnpeJ r £ (v £ ,(p £ ) < T(y) . 



a 



6.3 Proof of Corollary 1.2: 

Let (p E X with J-{(p) < +oo. From the upper bound inequality in Theorem 1.1, there is 
a sequence (wi, e ,U2,e) £ H such that 

limsup£(£ £ (wi i£ ,w 2 , £ ) - £ £ (%)) < ^(<y5) • 
Since (ui j£ ,u 2>£ ) minimize S £ in H, the previous inequality yields 

limsup£(£ £ (wi, e ,W2, £ ) - E £ (r) £ )) < F(ip) , (6.14) 

so in particular (ui, e ,ti2,e) satisfy (1.16). Hence, from the compactness and the lower 
bound inequality in Theorem 1.1, there is y? G X and a subsequence (iii, e ', W2,e') with 

liminf e' (£ £ >(u l£/ ,u 2 , £ ') - E £l (r) £ >)) > J"(</?) . 
This inequality is verified for every subsequence of (ui t£ ,u 2j£ ), so we have 

liminfe(5 e («i >e ,t42, e ) - E £ (rj £ )) > F{ip) . (6.15) 

From (6.14) and (6.15), we obtain 

F((p) > limsupe£ £ (Mi i£ ,M 2 , £ ) > liminf e£ £ (wi i£ , w 2 , £ ) > ^(y?) , (6.16) 

£ ^o £ ^° 

so J-(ip) = infx T. Taking (p = ip in (6.16) yields 



\\me(£ e (u^ e ,U2,e) - E £ (rj £ )) = inf T . 

e^O X 



u 
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6.4 Proof of Corollary 1.3 

We start proving that when a,\ is not to close to or 1, the minimizers of J 7 in X are not 
radially symmetric. We show that for any radially symmetric i/jfl, J 7 if) > J~((pd s ), 
where the support of (fids G X is a disk sector. We first prove this for functions such that 
{(/? = 0} is a disk or an annulus. Then, we generalize by induction the result to radial 
functions such that {</? = 0} is composed of a finite number of connected components. We 
conclude then by approximating any radially symmetric (p G X by this kind of functions. 

We recall that p is given in (1.3) and that X is the space of functions if G BVi oc (V; {0, n}) 
such that 

p — (X\. (6-17) 

{<p=0} 

If ifds G X is such that {ipd s = 0} is a disk sector, we easily compute 

8a ~ 16 ' 

For < R~ < R + < Awe denote A(R~ , R + ) the annulus of center the origin, inner radius 
R~ and outer radius R + . 

If ip a G X is such that {</? = 0} = A(0,R a ) and f A , QR \P = ct, then R a = A(l — \fa) x l 2 
and 

^ = /(«), (6-18) 

OCT 

where / : [0, 1] — >■ R + is the concave function f(a) = (1 — a)^ 4 (l — y/1 — a) 1 / 2 . We see 
that there exists 

S « 0.1486 
such that if a G [#o, 1 — So], then /(a) > 3 /i6. 

Proposition 6.1. If a± G [#o,l — So], then the minimizers of T in X are not radially 
symmetric. 

Proof: (Step 1) Let R G (0, A) and consider <p{ G X such that {ipf = 0} = A(0, i?). From 
(6.17), we have that J r (v?f)/8<r = /(«i) so (6.18) yields 

F(<pi) > Ffrds) . (6.19) 

Since a 2 = 1 — «i G [5 , 1 — 5 ], the similar inequality holds if {y?f = 0} = A(i?, A). 

Consider now <p\ G X such that {<£>" = 0} = A(i?i, i? 2 ), with < R\ < R 2 < A. Writing 



we compute 



Pi= P, p2= P and (3 3 = p. 

lA(0,Ri) JA(R 1 ,R 2 ) JA(R 2 ,X) 
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OCT 

From (6.17), we have that (3 2 = «i and /?i + /?3 = «2 so 

^^ = m) + m + «o . 

SO" 

The right hand size of the previous equality is a concave function of /3\ and the value of 
f3\ may vary between and a 2 . If f3\ = then ^{ip'D/Sa = f{ot\). If /?i = 0:2; since 
«i + a 2 = 1 we find J r (</?")/8(j = f{ot 2 ). We derive 

^(^) > F(ip ds ) . (6.20) 

(Step ,2) Let n G N* and consider <y9 n G X such that 

n 
{fn = 0} = |J % , 



with A 2 j = A(R 2 j,R£j) and 



< R 2 j-2 < -^2j-2 < ^2j < ^2j ^ ^ 



for 2 < j < n. We write /3 2j = J Aaj P, Pi = J A{0 ^) P, /Wi = L«,a) P and 

/?2j+l = / P 



'A(R 2j ,R 2j+2 ) 

for 1 < j < n — 1. Notice that we allow A(0, i?^) or A(i?^, A) to be empty but this only 
implies that (3i = or (3 2n+ i = 0. With this notation, we have 



n 


= «!, 




= a 2 



(6.21) 



and 



j=i \i=i / 



By induction, we are going to prove the following property: 



n 



(V„) V/3i, • • • , (3 2n+ i G [0, 1] such that ^ (3 2i = a-y and ^ # 

9n(l3l,--- ,02n) > 



2i+l — &2 , 
i=l i=\ 

^{fds) 



8a 
If n — 1 we are in one of the three cases analyzed in Step 1, so (6.19) and (6.20) yield (Pi). 

Let us assume that (V n ) holds and consider j3±, ■ ■ ■ , f3 2n +3 G [0, 1] such that 
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rc+1 ra+ 1 

^/3 2 i = «i and ^/3 2 j+i = «2 • (6.22) 

We have 

2n / j \ /2n+l \ /2n+2 

£n+l (/?!, " " " , /3 2 „ +2 ) =E/EA+/E^+/E^ 

j=i \i=i / V i=i / \ j=i 

The right hand side of the previous equality is a concave function of /?2n+2- The value of 
/3 2 n+2 may vary between and a.\. Suppose first that /?2n+2 = 0. Then, defining 

J3j = (3j if j = 1, • • • , 2n and /3 2 n+i = /Wi + /W3 , 

the /3j's satisfy (6.21) and we have 

2n / j \ 

g n +l(/3l,--- ,/?2n+2) > 5^/ ( ^2 Pi J =fi'n(/3l,--- , #2n) • 

Hence, (V n ) yields # n +iC#i, • • • , /W2) > F{(Pds)/%v. 

Suppose now that /?2n+2 = ot\. From (6.22) this implies 02j = for every j = 1, • • • , n. 
Then, defining 



2n+l 

°2n+3 , 



Pi = ^2 Pi > P" 2 = P 2n + 2 and ^ 3 = ^ : 



the fa's satisfy (6.21) and we have 



'2n+l \ /2n+2 



9n + l(Pu--- ,P2n + 2) > f [Y, & } + f [Y, & 

= f(Pl)+f(Pl+P2) 

= 9i(Pi,p2). 

Hence, (Pi) yields g n +i(Pi, • • • , P2n+2) > J~(<Pds)/8&- We derive that the result holds for 
all the possible values of /3 2 „ + 2- 

We have proved that if tp n e X is radial and its support has a finite number of connected 
components, then 

^(y? n ) > F{(p ds ) ■ (6.23) 

(Step 3) Suppose now that if G X is a radially symmetric function such that {ip = 0} has 
an infinite number of connected components. Since tp has locally finite perimeter in V, 
{ip = 0} is the union of a countable family of disjoints annuli. We write 

jez 
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with A 2 j = A(R 2 j, RJj) such that 

< i? 2 " < i2+ < R 2j+2 < R+ j+2 < X . (6.24) 

For every n G N, we define a function <p n : V — > {0, n} by 

{<^ n = 0} = (J A 2i (J i 2n+2 (J i_ 2n _ 2 , 

j=-n 

such that 

A 2n+2 = A(L+, A) and i_ 2n _ 2 = A(0, L~) 
with L~, Lt\ to be chosen next. If (L~, L+) = (0, A), then (6.24) gives 

p=J2[p < [ p- 

Similarly if (L~,L+) = (RZ 2n , #£J, then 

p=£/ p+J2 [ + p+H f + p> f p- 

Hence, by continuity there is a pair (L~, L+) G (0, RZ 2n ) x (ifcjrn A) such that A- =Q , p = 
// =oi'° = ai- Clearly <p n G BVi oc (V), so <p n G X. Moreover, (6.24) yields 

lim L~ = and lim L+ = A . (6.25) 

We have 



^) = W p 3/2 ^ , 



and since p is radially symmetric 

n » 

•Ffan) = E/ p 3/2 ^ 1 + 27r(p 3 / 2 (L+)L+ + p 3 / 2 (L;)L;). 

j= _ n JdB(0,R+) 

From (6.25), the last term in the previous equality goes to zero as n — > +oo, so lin^oo J-(ip r 
F{(p). Hence, since {ip n = 0} has a finite number of connected components, (6.23) yields 
J- {(f) > J-{<4>ds)i which ends the proof. 

□ 

Proof of Corollary 1.3: Suppose that a.\ G [So, 1 — So] and that {(ui j£ ,u 2:£ )} £>0 is a 
sequence of radially symmetric pairs such that (ui j£ ,u 2t£ ) minimizes S £ under the mass 
constraints (1.5). Then, tp E defined by (3.7) is also radially symmetric. Consider <p £ fi, the 
restriction of (p £ to a slice of T> passing through 0. From Proposition 5.1, <p £) o belongs to 
SBVioc([0,X] ; {0, 7r}) and converges in L} oc ([0, A]) to <po- Hence, ip £ converges in Lj oc (V) 
to the radial function tp given by tp(x) = (po(\x\). From Corollary 1.2, we know that (p 
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minimizes T over X, which yields a contradiction with Proposition 6.1. 

□ 



7 Appendix 

We end this article given the proof of Lemma 4.2, which is essentially the same of Lemma 
4.3 in [10], which in turn is a generalization of Lemma 1 in [26]. For completeness we give 
here the details of the proof. 

Proof of Lemma 4.2: (Step 1) Suppose first that T>C\A and T>\A have both non empty 
interior and let 



B(x 1 ,5)cVnA and B(x 2 ,5)cV\A. (7.1) 

We first approximate A by sets of finite perimeter in D. For k > 2 we define T> k = 
V n 5(0, A(l - l /k)) and A' k = An V k . We have that d*A' k C (d*A n V k ) U dV k , so 

3/2 ^< / p^d^H- / o i/2 dV} 



p 1 " an s / p 1 an + / p 



Using Lebesgue dominated convergence theorem, the first term in the right hand side of 
the inequality converges to f d A p^ 2 ctH 1 < +oo. The definition of V k and (2.8) yield 



f /2A \ 3 / 2 

\ m P ,2 dU l < \\ P \t[l {&Dk) V}{&D k ) < (— ) U\dV) = 0,^(1) . 

Hence, 

lim / p^dH 1 < [ p i/2 dV}. (7.2) 



k ^ oc Jd t A' k Jd,A 



(Step 2) Since A' k has finite perimeter in V, it can be approximated (see the proof of 
Lemma 1 in [26]) by open bounded sets A k , such that 

C\A k AA' k ) < I (7.3) 

A' fc C4 + 5(0, i/fc) and A k C A' fc + 5(0, i/ fc ) (7.4) 

^(SAfcnSP) = 0. (7.5) 

The definition of A' k and (7.3) imply (i). Using (7.1) and (7.4), for large enough k we get 

B(xi,6)cA k and B(x 2 ,5) C V\A k . (7.6) 

Moreover, using (ii) from Proposition 2.3 in [10] and the fact that A k belongs to a sequence 
A k such that ||-4fc||Bv(x>) - > \\A k \\BV(D) as n — > 0, we have 

P^\Dl Ak \< f p^\Dl A ,\ + \. (7.7) 
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Also, the definition of A' k and (7.3) yield 



Ik := / P~ / P = Ofc-K»(l) ■ (7. 

'A fc ./A 



(Step 3) Now, we set 



A k \B(x 1 ,r ltk ) if 7fe > 

.4/, = { _ A k if 7fe = o 

AfcUS(xi,r 2lfc ) if 7fc<0 
where ri^ and r 2i fc are chosen to satisfy 



P = P = lk- 

B(zi,n,fc) J B(x2,r 2 ,k) 



Since r (->• J* B , s p is continuous and decreasing for r G (0,5), r lfe and r 2ife are unique 
and tend to zero as k — > oo. Then, we derive from (7.6) and (7.8), for large enough k, 
that 



P = / P~lk= I p. 

VnA k J A k J A 



Moreover, from (7.5) and (7.6), we have 'H 1 (dAk) = for k large enough, so (ii) is proved. 
Using again (7.6) we obtain 

/ p 3/2 \Dl Ak \ < f p 3 / 2 \Dl Ak \ + \\p\\ 00 H 1 (dB(x 1 ,r 1:k )UdB(x 2 ,r 1:2 )). (7.9) 
Jv Jv 

Hence, using (7.7), we obtain 

/ p /2 \Dl Ak \ < [ p /2 \Dl A , k \ + cw,(l) , 
Jv Jv 

so (7.2) gives 

limsup f p 3 ' 2 \Dl Ak \< f p 3/2 \Dl A \ 

k-^oo Jv Jv 

We have proved (iii). 

(Step 4) We now remove the condition that T>(lA and T>\A have no empty interior. First, 
we notice that £ 2 (V fl A) = and C 2 (V\A) = are not possible because of the mass 
constraints in (1.8). Hence, there exists Xi a point of density of T> D A and x 2 a point of 
density of T>\A. Consider the function 

$(<M 2 ) = / p- p, 

J A 12 J A 

where A\ 2 = A U B(x±, Si)\B(x 2 , S 2 ). Since p > in V, for any 5 > we have 

$(<5,0) >0 and $(0,5) <0. 

Since $ is continuous, there is t — tg G (0, 1) such that $(£<5, (1 — t)S) = 0. Define 
As = A U B(xi, (1 — t)S)\B(x 2 , tS) and <p = ttI^. Both DnAj and T>\A$ have no empty 
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interior and f A p — f A p. Moreover C 2 (AgAA) — > as 5 — > 0, and using an inequality 
similar to (7.9), we get 



limsup f p 3/2 \Dl Ag \ < [ p 3/2 \Dl A \ . 

8^0 Jv JV 



Finally, for each Ag we apply the construction from steps 1-3 and conclude thanks a di- 
agonal argument, see Corollary 1.16 in [7]. 

□ 
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